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Abstract 

In this paper we derive lower bounds in minimax sense for estimation of the in- 
stantaneous volatiUty if the diffusion type part cannot be observed directly but under 
some additional Gaussian noise. Three different models are considered. Our technique 
is based on a general inequality for Kullback-Leibler divergence of multivariate normal 
random variables and spectral analysis of the processes. The derived lower bounds are 
indeed optimal. Upper bounds can be found in [18) . Our major finding is that the 
Gaussian microstructure noise introduces an additional degree of ill-posedness for each 
model, respectively. 
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1 Introduction and Discussion 

Let Xt = a {s) dWs^ where (Wi)^^^^ denotes here and in the following a standard 
Brownian motion. Consider the model 

yi,n = Xi/n + rei^n, i = l,...,n, (1.1) 

where e M (0, 1) , i.i.d. Throughout the paper (Vl't)jgjQ and (ci^nj • • • i ^n,n) sre as- 
sumed to be independent. Further a is an unknown, positive and deterministic function 
and T > is a known constant. 
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In the financial econometrics literature variations of model (jl.ip are often denoted as 
high-frequency models, since (VFt)^g[Q is sampled on time points t = i/n corresponding to 
short time intervals 1/n which can be of the magnitude of seconds, nowadays. There is a 
vast amount of literature on volatility estimation in high-frequency models with additional 
microstructure noise term (see Barndorff-Nielsen et al. [1] , Jacod et al. [15] , Zhang [23] and 
Zhang et al. [24] among many others) . These kinds of models have attained a lot of attention 
recently, since the usual quadratic variation techniques for estimation of a^{s)ds lead to 
inconsistent estimators (cf. Zhang |23j). Brown et al. [[B] studied low- and high-frequency 
volatility models by means of asymptotic equivalence. Recently, Reifi |20) has shown that 
model is asymptotically equivalent to a Gaussian shift experiment. 
Closely related to model (jl.ip is 

Yi,n = a {^) Wi/^ + rei^n, i = l,...,n. (1.2) 

This model can be regarded as a nonparametric extension of the model with constant a, t 
as discussed by Gloter and Jacod [11], [12] and for variogram estimation by Stein [2lj. As 
a further natural modification of (II. ip . we consider 

Yi^n= / Xgds + rei^n, i = l,...,n. (1.3) 
Jo 

For constant a, the process 

( f Xsds) = ( [\t - s) a (s) dWs) 

is called integrated Brownian. In this case, model (jl.3p has been used as a prior model for 
nonparametric regression (see e.g. Cox [8]). 

All models have a common structure: They might be interpreted as observations coming 
from a particular Volterra type stochastic integral, i.e. Jq K (s, t) a (s, t) dWg under addi- 
tional measurement noise. In this paper we derive lower bounds for the models (jl.ip - (jl.3p . 
We stress that the treatment for general K is an interesting but rather difficult task. 
More precisely, we derive minimax lower bounds for estimation of the instantaneous volatil- 
ity, i.e. o"^ as a function of time, with respect to -L2-I0SS. One of the key steps is to consider 
the estimation problem after taking finite differences which is typical for variance estimation 
(see Brown and Levine [1] or Munk et al. |T6]). Usually in nonparametric regression, lower 
bounds are obtained under an independence assumption on the observations. In order to 
deal with dependend data, we introduce a new bound of the Kullback-Leibler divergence 
which might be of interest by its own. The lower bounds then follow from a standard 
multiple testing argument together with this bound and the control of the eigenvalues of 
the covariance operator. 
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In nonparametric variance estimation, n""/'-^""''^) is well-known to be the minimax rate of 
convergence given Holder-smoothness a (for a definition see ()2.2|) ) of the variance and a 
sufficiently smooth regression function (see Brown and Levine [4J or Munk and Ruymgaart 
|17j). We show that for the microstructure noise models (jl.ip and ()1.2p the lower bounds are 
^-a/(4o+2) ^^^^ £qj, ]-Qodel ()1.3p . 77,-"/(8°+4)_ jf (j exceeds some minimal required smoothness 
these rates are shown to be also upper bounds for models (jl.ip and (|1.2p (Munk and 
Schmidt-Hieber [18j). 

For constant o", i.e. Yi^n = (^Wi/n + 8Ta^n~^^^ is the optimal asymptotic variance 

for estimation of a^. First this has been shown in a more general setting by Gloter and 
Jacod but can be also derived as in Cai et al. [7] with respect to minimax risk via 
the information inequality method (see Brown and Farrell [3] and Brown and Low [5j). 
In order to explain the optimal rate n"^/^ it is tempting to think that this is related to 
the pathwise smoothness of the Brownian motion, since optimal estimation of a Holder 
continuous function with index 1/2 results in a n~^/'^ rate of convergence. However, this 
reasoning is in general not true since the smoothness of integrated Brownian motion is 
arbitrarily close to 3/2. By the same argument we should obtain in model (|1.3p an n~^^^ 
rate if a constant. This contradicts the obtained lower bound n^^/^. 

Indeed, we find it more properly to look at these models from the viewpoint of statistical 
inverse problems. The eigenvalues of the Brownian motion and integrated Brownian motion 
covariance operators behave like 1/i^ and 1/i^, respectively (see Freedman [9]). So the max- 
imal frequency i for estimation of a is reached when l/i^ ~ 1/n (for Brownian motion) and 
~ 1/n (for integrated Brownian motion), i.e. in the first case we may use O (ni/2) fre- 
quencies and in the second case O (n^/^) resulting in the reduction of the rate of convergence 
by a factor of 1/2 and 1/4, respectively. Motivated by this heuristics, we conjecture that the 
optimal rate of convergence for the kernel K (s, t) = {t — sY ,q £ [0, oo) and a (s, t) = a (s) 
is n~"/[('^'^+'^^^'^""'"^^l. This implies that the rate of convergence decreases as the order of the 
zero of K increases, or equivalently, as smoother the path of K (s, t) a (s) dWg becomes. 
For a further result in this direction (g G (0,1/2), a constant) see Gloter and Hoffmann 
[To]. Note that for q € [1/2, 1) U (l,oo) the spectral decomposition of the integrated Brow- 
nian motion is not known and hence our strategy of proof cannot be applied. More general 
techniques are required. 

Note finally that the lower bounds still hold if we consider generalizations of the models. 
For instance consider model (jl.ip and allow a to be random itself. Assume further that 
e = (ei,n5 • • • ,era,n) is general microstructure noise, i.e. white noise with bounded moments 
(see Huang et al. [14J or Zhou [25]). In this model the lower bound for estimation of 
the instantaneous volatility is still of course. To show that this is indeed the 

upper bound is current research by Marc Hoffmann and the authors, where the estimator 
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is constructed by wavelet techniques and methods as in Brown et al. [2] are employed. 



2 Results 

The next Lemma might be of interest by its own and can be applied to various problems 
in variance estimation where we do not have independent observations. The Lemma can 
be viewed as a generalization of Lemma 2.1 in Golubev et al. |13) for matrices with non- 
uniformly bounded eigenvalues. See also inequality (3.8) of Reifi [19]. For our purpose it 
is required to allow eigenvalue sequences tending to and oo. The proof of our minimax 
results follows standard arguments (see Tsybakov [22]). Note, however, that the under- 
lying dependency structure of the process causes severe technical difficulties and requires 
very sharp bounds for the Kullback-Leibler distance in these models as given in the next 
lemma. Recall that the Kullback-Leibler divergence between two probability measures P, Q 
is defined as 



(P. <?):=/ log QdP. 



whenever P <^ Q and -|-oo otherwise. 

Lemma 2.1. Let X ~ A/" (/i, Sq) and Y ^ M {fj,, Si) be n-variate normal r.v's with expecta- 
tion fj, and covariance Sq and Si, respectively and denote by Px and Py the corresponding 
probability measures. Assume < CSq < Si for some constant < C < 1. Then it holds 
for the Kullback-Leibler divergence 



dK(PY,Px) < 



4C2 



— 1/2 —1/2^ ||2 

Sg (Si-So)Sq ^ < ||So Si - , (2.1) 



where /„ denotes the n x n dimensional identity matrix and \\-\\p is the Frobenius norm, 
i.e. for a matrix A, \\A\\p := tr^/'^ (vlvl*). 

Proof. Note that 

dK{PY,Px) = \ (^log (^det (s~'/'SiS~'/')"'^ +tr (s-'/'SiSo'/') - . 

Introduce S := Sq ^^^SiSq and note that S is positive definite. Further, by assumption 
< CSq < Si and hence Sq ""^^^SiSq ^^"^ > Cln- This gives Wi := A, (S) > C, i = 1, . . . ,n. 
Recall det (S) = HiLi ^« (^) ~ Z]r=i '^i- Hence 

^ / n n \ 1 " 

dK{PY,Px) = 2 - (^^) + ^Wi-n] =-J2iwi-log{l + Wi-l)-l). 

\ 1=1 1=1 / i=l 
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Assume x > C — 1. Expand — log (1 + x) = —x + ^2 (1 + ^) j x^ for a suitable |^| < 
For C-l<x<Owe have - log (1 + x) < -x + x^/ (2C^) and for x > 0, - log (1 + x) < 
—X + x^/2. Therefore by Lemma lA. 41 



1 " 

dK{PY,Px)<^^{wi-l 



1=1 



4C72 



.-1/2 



(Si 



So) S-^/' 



□ 



Remark 2.1. T/ie assumption < CSq < Si can be relaxed at the cost of an additional 
symmetrization term in \2.1\) . More precisely, assume in Lemma \2. 1\ instead ofO< CSq < 
Si that < So, Si holds. Then we have 

dK{PY,Px) <^\\^^'^l-ln\\l + \\\^i'^0-ln\\l. 

Next we prove a lower bound for Holder continuous functions cx^ and r^. For this 
we need some notation. We write [x] := max^g^jz < x}, x G M, the integer part of x. 
Let < / < n < cxD some constants. The class of uniformly bounded Holder continuous 
functions of index a on the interval I is defined by 



C'' {a,L) := {a,L, [l,u]) := |/ : /(^^ exists for p = [a] , 



(2.2) 



< L\x 



\a—p 



y\ ^ , Vx, y G I, 0</</<u< cxd| . 

For any function g we introduce the forward difference operator Aig := g {(i + 1) /n) — 
g (i/n) . log() is defined to be the binary logarithm and we write Mp and Bp for the space of 
px p matrices and px p diagonal matrices over R, respectively. The i-th largest eigenvalue 
of a Hermitian matrix M is defined as Aj (M) . 

Theorem 2.1. Assume model and a > 1/2 or model il.^) . a > 1. Then there exists 
aC > (depending only on a,L,l,u), such that 

liminfinf sup E ( lla^ - fj^ll^ ) > C. 



Proof. Besides working with the observations li 1^ ^ directly, we consider the suffi- 



cient statistics 



Ay 
Ay 



{y^l,ni ^2,n ^l,ni • • • , "^n,n "^n—l,n) , 
"^l.ni ^2,n ^l,ni • • • i ^n,n ^n—l,n I • 



Let AiW := VF(j+i)/„ - and Ajei,„ = ei+i,„ - ej,n- We obtain for i = 1, . . . 



n 



Yi. 



Yi—l^n 
Yi—ln 



a 



i/n 
{i-l)/n 

i 
n 



a{s)dWs + rAi_iei_i,„, 
Ai_il^ + (A,_io-) VF(i_i)/„ + rAi_iei_i,„, 



(2.3) 
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where eo,n := 0, Yo.n := and lo,n := 0. We may write Ay = Xi + and Al" = 
X[ + X2 + Ri, where Xi,X[,X2 have components fl^.^-^y^a{s)dWs, a{i/n)Ai^iW and 
rAi_iei_i,„, respectively and = (Ai_ic7) 

Let us first consider model p.2p . a > 1. By the ease of brevity, we simply point out in a 
second step the differences to model (jl.ip . The idea is to prove the lower bound by a multiple 
testing argument. Explicitly, we apply Theorem 2.5 in Tsybakov [22) . The construction 
of hypothesis is similar to the one given in [22], Section 2.6.1. We write o"mirnCmax for the 
lower and upper bound of o"^, respectively, i.e. ci^ G {ci,L, [crmim (^max])- Without loss of 
generality, we may assume that Cmin — 

1. Let 

K:M^M+, K{u) = aexp(— L-^V||2„|<i}, (2.4) 

V 1 - (2n) / 

where a is such that K £ C {a, 1/2) . Further for some c > 0, specified later on, let m := 

[2-icni/(4a+2) + !],/!„ = {2m)-\ tk = K{k- 1/2) + 1/4, 

<Afc (t) := L/i^i^ (^^^) , A: = l,...,m, tG[0,l]. 

Define 17 := {w = (a;i, . . . , w^) , G {0, 1}} and consider the set 

£ ■= [al (t) : (t) = 1 + YlT=i^kMi)^^ e ^} • Then it holds Vcj,cj' G f7 

IkS - -I'Wl = f ('^'(^) - ^l'it)fdt = L^e+i ||K||2 p (^,^') , (2.5) 

JO 

where p{u),u)') = X^^Li | is the Hamming distance. By the Varshamov-Gilbert 

bound (cf. Tsybakov ^2]) there exists for all m > 8 a subset {loq, . . . ,ujm} of 0, such that 
^0 = (0, . . . , 0), p {uji, LOj) > m/8, V < j < < M and M > 2"*/^. Define the hypothesis 
Hi for i = 0, . . . , M by the probability measure Pi induced by af^ := cr^. By Theorem 
2.5 in Tsybakov [22] the proof is finished once we have established 

(i) cTj- „ E C (a, L) , 1 < < cxmax 

(u) ||ct2„ - a]J\^ >2s> cn-°/(4-+2)^ i / j, c> 0, (2.6) 
1 

(m) — X^^i^l-Pj'^o) <KlogM, j = l,...,M, K<1/10. 



M 



(i) is obviously fulfilled for sufficiently large n. By (12. 5p it follows for m > 8, 
16~^L/i^ ||i^||2 and hence (zi). (izi) We apply Lemma [2. II in combination with Lemma 
Let Ilfc G D„, k = 1, . . . , M with entries (H^) ■ = (Tfc^„ {i/n) dij. For the observation vector 



> 

_2 

X2I 
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Y = 

( ^i,n) • • • 1 Yn^n ) j have Y ~ J\f (0, S^) under Hj^, k — 0, . . . , M , where 



iAj 



n 



Hi 



i,j=l,...,n 
i A j~ 



n 



1,...,M. 



i,j=l,...,n 



Because of cta;,^ > 1 it fohows \ak,n {x) - (yk,n {y)\ < crl^n (^) ~ ^fc,n iu) ^ k - y| , i-e. 
a £ C {1, L) . By Lemma lA.21 < (2 + 12L^) ^ Sq < S^. These inequahties remain vahd 
under any invertible hnear transformations of Y. Hence if we denote by the covariance 
of Ay under Hk {k = 0, . . . , M) it fohows < (2 + 12L'^)~^ Sq < and we may apply 
Lemma O with C = (2 + 121"^)'^. Hence 



(2 + 12L2)2 _^ 2 
dK{Pk,Po) < ||Sq T,k — In\\p k = I, 



M. 



Let for 1 < i, j < n, the matrix ^ be defined by 



2 for i = j and i > 1 
— 1 for \i — j\ =1 
1 for i = j = 1 
else 



(2.7) 



and let := H^. - /„, k = 1,...,M. Clearly, Tk < Lh^\\K\\^. We abbreviate the 
covariance of two column vectors X and Y as the matrix with covariances of XY^. Then 
we have the explicit representations 

1 



^0 



-In + tM, 



n 



1 



So + -Tk + -Tl + CovH, {X[, + Covj^, {Ri,X[) + Cov^, [Ri] 



n 



n 



where the subscript Hk means that these covariances are taken with respect to the prob- 
ability measure induced by Hk- We remark that due to Sq > In/n it holds Ai (Sq ^) = 
A^-*^ (So) < A^"*^ {In/n) = n. This yields using Lemma [A. 31 (i) 



dK iPk,Po) < 5-&±^^ (± W^o'^kWl + ^ W^o'^lWl + ||So ' CovH, {X[,R, 



|2 



+ \\^,'CovhARuX[)\\1 + \\j:^'CovhARi) 



<5il + 6L')'((4L'\\K\\l ' 



AL'\\K\\l^ + L^\\K\\lhl^]hl^n-^\'^^ 



+ r? \\CoYH^ {X[,Ri)\\l + n' ||Covh, X() ||^ + ||Covj^, (i?i)||^) . 
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The remaining part of the proof is concerned with bounding these terms. We make use of 
the properties on Frobenius norms collected in Lemma lA.41 and obtain 

2 



Let E £Mn given by 



(iAi)-l 



n 



i,j=l,...,n 



1 if j > i 
otherwise 



Also let Allfc G IDn, /c = 1, . . . , M with entries (Allfc)^ j = (Aj_icrfc^„) 5ij, where o"fc^„(0) = 0. 
Then Covh^ {X[,Ri) = n~^IlkE (Allfc) and for n large enough 

llCovi^, {X[,Ri) \\l < 4n"2 tr (^E (AUk)'^ < AL'^n''^ \\E\\% < AL'^n'^. 

With the same arguments we can bound ||Sq ""^ Cov/^j, Further we have for 

j = 1, . . . , M using Lemma lA.ll 



i=l 



for a constant Cr only depending on r. This gives 
1 *^ 

< 5 (1 + 6L^)' ((4L2 lli^ll^ + \\K\t C) n'/^hl^Cr + + 8L^ 

< Cr,L,\\K\\^n^''^h^n < Cr,L,||i^||^c"^"~^2m < K log M, 
where C^,L,||_ft-||^ is independent of n and the last inequality holds if 

/ _ Nl/(2a+l) 

C > \ lQCr,L,\\K\\^l^ j , 

using the Varshamov-Gilbert bound. 

The proof for model (jl.ip and a > 1/2 is almost the same as for Theorem 12.11 So we 
only sketch it here. Note that Lemma 12.11 can be applied directly without use of Lemma 
IA.2I The construction of hypothesis is the same. Let 11^ G D„, k = 1, . . . , M defined by 
(Ilfc)^ ^ := fl^-^iy^ cr'i n i^) '^^ ^^"^ '■= — In- With the notation as in Theorem 12.11 the 
problem can be reduced to a testing problem where we have to test the hypothesis that 
a centered random vector has covariance matrix Yiq = In + A against the M alternative 
covariance matrices Sfc = Sq + Ffc, k = 1,...,M. With maxi<j<„ maxi<fc<jv/ (rfc)j ^ = 
O {h") and Lemma [XT] the Theorem follows. □ 
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The proof of the lower bound for estimation of cr^ in model (jl.3p differs from the 
previous one. Instead of using first order differences we transform the data by taking 
second order differences. The key step is to show that for constant a this is "close" to the 
model where we observe = crn~^/'^rji^n + T(,i^n, i = 1, ■ ■ ■ ,n, iji^n ~ A/" (0, 1) , i.i.d. Here, 
= iCi,n, • • • ) Cn,n) is & particular MA(2) — process, independent of ?? = (r/i,„, . . . , r/„,„) . 

Theorem 2.2. Assume model il.3\) and a > 1/2. Then there exists a C > (depending 
only on a,L,l,u), such that 

lim inf inf sup 

Furthermore, there exists aC > (depending on cTmin; o'max/' such that for constant a and 

< 0"min ^ O^max ^ OO 

lim inf inf sup E (n^/^ (a"^ - a^) ^) > C. 

n— >-oo 5-2 2^r i V / 

Proof. Except for changing the definition of m := ^2~^cn^^^^"~^'^^ + l] to m := [2~^cn^/(®°'^^^ + l] 
we construct the same hypothesis as in the proof of Theorem 12.11 In order to prove the first 
part of the statement it remains to show (iii) in (|2.6p . We consider second order differences, 
i.e. 

Y* := AtiF,, Fo := 0, i = 2, 3, . . . , Y{ := V2Yi, 
where A? = Aj o Aj. Note that for i > 2 

pi/n / • \ Ai-\)ln f j-2\ 

Y*= [--s)a{s)dWs+ is ]a{s)dWs + TA^ei 

J{i~l)/n \n J J{i-2)/n V n J 

Let Y* := {Y^ , . . . ,Y*) . Obviously, this is equivalent to observing Y . The covariance of Y* 

under hypothesis Hk is denoted by T,k,k = 0, ...,M. Since by construction o"q„ < 

it follows from elementary computations that — Sq is the covariance of Y* under a = 
\ 1/2 

0"^ ^ — fJo „ ) and r = 0. From this we conclude that 

Sfc - So > 0, (2.8) 

i.e. So < Sfc. Note that Var {X + Y) < 2Var X + 2Vary (in the sense of Loewner ordering, 
see also Lemma |A.3|, (iii)). This yields together with ()2.8p 

Sfc — So < r, 

where F is diagonal with entries 

^'h,j- 3^3 "m' 
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and Sij denotes the Kronecker delta. Let denote the probabihty measure of under 
iffc. Due to So < Sfc we may apply Lemma |2. II and obtain 



1 



-1/2 ^ ^ V-V2 

4 

Direct computations give 



dK{Pk,Po)<^ ^o ' i^k-^o)^o '' ^< — II VII;- (2-9) 



2 ^ WL^hl-\\K\\l II ,,,2 



n 

where A is as defined in (j2.7p 



1 if i = 1, j = 2 or i = 2, j = 1, 
otherwise. 



and V2 is symmetric and {V2)ij 7^ only if i,j < 3. Obviously, the smallest eigenvalue of 
Vi is —1. Hence we can estimate by Lemma lA.lt 

So > ^In + {A^ + V2) . 

Since V2 has only non-zero entries in the first three rows and columns, it has only three 
non-zero eigenvalues. By standard bounds on eigenvalues (see Lemma [A. 31 (ii)) this allows 
to estimate for i > 3 

\n-i {A^ + V2) > Xn-i+3 {A^) + A„_3 {V2) = A^.+S (^) • 

Let Vn '■= [ra^/^] . Then for sufficiently large n by Lemma lA.ll 

IIVll^ = E^HV)<E^n%(^M+-^ E K'-.^ii^'+y^) 

2=1 i=l ^ ^ «=r„+l 



n 



« = »'n + l i=T-„ + l ' 



< 36n25/4 + J_ y 1 

„l/4 /. , 1 



Because the last part converges as a Riemann sum to a finite integral, we may find a 
constant Cr depending only on r such that ||Sq ^||^ < Crn^^/^. This gives with (12. 9p . 

1 *^ Ifi 2* 

-J2dK{Pk,Po) < -CrL'\\K\\lhl"n'/' < -CrL'\\K\\lc-"^-'m<KlogM 

k=l 

/ o \ 1/{2q+1) 

whenever c > (2^/9CrL^ W^Wlo '^'^ 



In order to prove the second statement of the theorem we consider two hypothesis ctq := Cmin 
and af := (Tmin + cn~^/^ and apply Theorem 2.5 in Tsybakov |22] for M = 2. Using the 
bounds from the first part, the remaining part of the proof is straightforward and thus 
omitted. □ 
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Appendix A Technical tools 

Lemma A.l. Let ^, Q"^ be as in (2.1) and (HHP; respectively. Then 



A„_.+i {A) = A„_.+i (Q-i) = 4sin2 ( i^LJll^^ > ^' 



Proof. The i-th eigenvector of Q^"*^ is given by Vi = (sin (xj) , sin (2xi) , . . . , sin (nxj)) , 
where Xi := {2i — 1) vr/ (2n + 1) , i = 1, . . . , n. The corresponding eigenvalues are 

'{ 2i - l)7 r' 
4n + 2 

If f j = . . . , t'i,n) is an eigenvector of so is Vi = {vi^n, ■ ■ ■ , an eigenvector of A 
with the same eigenvalue. Using xtt/2 < sin (xvr) whenever x G [0, 1/2] , we obtain 



Xn~i+1 {Q-') = 4 



sm 



4sin^ {xiTT/2) > x-Tvyi > 



4(2n + l)2 4n2' 



□ 



Lemma A. 2. Let cr > 1 be a Holder C(l, L) function and let E G IDn 6e a diagonal matrix 
with T,i^i = a (i/n). Further introduce Q := {i A j)^ T/ien 

(2 + 12^2)"^ Q < SQS 

in i/ie sense of partial Loewner ordering of symmetric matrices. 

Proof. Obviously 

2 for i = j and i < n, 



Let O e M„ given by 



= < 


-1 


for \i — j\ = 1, 




1 


for i = j = n, 







else. 


















1 for i = j, 




— 1 for i = j — 






else. 



(A.l) 



Note that = OOK We have with 



{Aia'^Y for f=i-l, 
(Ajo--!)^ for i = i + l, 
else 
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that 



2 ^ ^ 2 2 



As can be seen by direct calculations 



-Ai(T" 



Define M G 



by 







\ 




2A„_i(7-2 






-A^_lCT-2 




1 


for i = j, 




Aia-2 - 1 


for i = j 


-1, 





else. 





Due to MM* > and ^EM it follows 



-2\ 



I (AiCT 



(A„_iCT- 



2\2 



V 

where we used in the last step that 
1 1 



4L2 



n 



2 



0/ 



(t2 (x) 



^2 (y) 



(t(x) -cr(y)| 



1 



+ 



1 



<2L\x-y\, 



(A.2) 



(A.3) 



0-2 (x) 0- (y) cr (x) 0-2 (y) 

i.e. £7 G C (1, and o" > 1 implies (7^2 ^ (i^ 2L). 

Next we bound the eigenvalues of the symmetric matrix S by showing that the cor- 
responding characteristic polynomial Xs(^) does not have any zero in [2a;,cxD), where w = 



max, 



In order to see this introduce the notation := ( S ) juj and note for 
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t > 2uj 



I tU}~^ -s{l) 
-sil) ■•• 



V 



> 



— s (n — 1) 
—s (n — 1) tu]^^ 

I l + s{lf -s{l) 



-s{l) 



1 1 -^(1) 



v 



1 + s (n — 1)^ —s{n — l) 
-s{n-lf 1 

\ / 1 



-s {n — 1) 
1 



V 



-s (n — 1) 
1 



> 



and therefore Xs(^) > for t > 2uj. Because oi w < I? jv? this shows that 
From Lemma I A. II fohows l/(4n^) < A^^^ (Q) = A„ (Q^^) and hence 

L2 



(A.4) 



This gives with (IA.4P finally 



□ 



In the next lemma we collect some important facts about positive semidefinite and 
Hermitian matrices. 

Lemma A. 3. (i) Let A,B£ Win ore positive semidefinite matrices. Denote by Xi{A) 
the largest eigenvalue of A. Then tr{AB) < Ai(A)tr(i?). 

(a) Let A,B £ Win CLfe Hermitian. Then 

Xn-r-s{A + B) > Xn-r{A) + Xn-s{B). 

(Hi) LetA,B are matrices of the same size. Then 

A^B + B^A < A^A + B^B. 
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Lemma A.4 (Probenius norm). Let A G M„. Then 

(i) 

n n 

\\A\\l:=iT{AA')=Y,K{AA') = Y,'^h 

i=l i,j=i 

and whenever A = A^ also = (A), 

(a) It holds 

4tr [A^] < ||^ + < 4p||^. 

(in) Let A, B he positive semidefinite matrices of the same size and < A < B. Further 
let X be another matrix of the same size. Then 

\\X^AX\\p < \\X^BX\\p. 

Proof, (i) and (ii) is well known and omitted, (iii) By assumption it holds < X^AX < 
X*BX. Hence Af [X^AX) < Af (X^BX) and the result follows. □ 
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